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ABSTRACT 



We investigate the quality factor and RMS amplitude of the lower kHz QPOs 
from XTE J1701-462 a unique X-ray source which was observed in both the 
so-called Z and atoll states. Correcting for the frequency drift of the QPO, 
we show that, as in all sources for which such a correction can be applied, the 
quality factor and RMS amplitude drops sharply above above a critical frequency 
For XTE J1701-462 this frequency is estimated to be ~ 800 Hz, where the 
quality factor reaches a maximum of ~ 200 (e.g. a value consistent with the 
one observed from more classical systems, such as 4U 1636-536). Such a drop 
has been interpreted as the signature of the innermost stable circular orbit, and 
that interpretation is consistent with the observations we report here. The kHz 
QPOs in the Z state are much less coherent and lower amplitude than they are 
in the atoll state. We argue that the change of the QPO properties between the 
two source states is related to the change of th e scale heigh t of th e accretion 
disk; a prediction of the toy model proposed by iBarret et al.l ( 120071 ). As a by- 
product of our analysis, we also increased the significance of the upper kHz QPO 
detected in the atoll phase up to 4.8a (single trial significance), and show that 
the frequency separation (266.5 ± 13.1 Hz) is comparable with the one measured 
from simultaneous twin QPOs the Z phase. 



Subject headings: accretion, accretion disks, stars: individual: XTE J1701-462, 
stars: neutron, X-rays: binaries, X-rays: individual: XTE J1701-462, X-rays: 
stars 
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Introduction 



Kilohertz quasi-periodic brightness variations (kHz QPOs) have been detected from 
more than twenty low-mass X-ray binaries containing a weakly magnetized neutron star 



( Ivan der Klia 120061 ). A pair of QPOs is often detected, with a frequency separation that 
typically changes by tens of Hertz as the individual QPO frequencies change by hundreds 
of Hertz. In lower luminosity systems (often called atoll sources after the tracks they make 
in color-color diagrams), the lower frequency QPO can reach quality factors (Q= v/Av) 



Barret et al. 


2005a 


b. 


2006) 



both the lower and up per kHz QPOs are broader features, with maximum Q of a few tens 
(IBoutelier et al.l |2010| ). In lower luminosity sources, for which the quality factor of the 
lower kHz QPO can be measured over its entire frequency span after correcting for the 
frequency drifts, Q increases with frequency until it reaches a maximum around 800-900 
Hz, beyond which a sharp drop-off is observed. The presence of the drop in many systems, 
and its reproducibi lity in a given system independent of count rate and spectral hardness 



(IBarret et al.l 120071 ). led to the suggestion that it may be related to the exi stence of an 



inner most stable circular orbit, a key prediction of strong field general relativity ( IBarret et al. 
2006h . 



XTE J 170 1-462 is a unique X-ray source, which first behaved like a Z source at high 
luminosity, and which later behaved like an atoll source at much lower luminosities. The 
source was clo sely monitored with the RXTE Proportional Counter A rray during its 2006- 



2007 outburst (IHoman et a 



20071 : iLin et al.ll2009bl ; IHoman et al.ll2010h . Type I X-ray bursts 



were studied by lLin et al.l ( 12009a| ). The three bursts observed occurred as the source was in 
transition from the typical Z-source behavior to the typical atoll-source behavior, at ~ 10% 
of the Eddington luminosity. N o significant burst oscillations in the range 30-4000 Hz were 
found during these three bursts. ILin et al.l (|2009al ) also derive a distance estimate of 8.8 ±1.3 
kpc from two radius expansion bursts (the latter two of the three). In both states, kHz QPOs 
have been reported, with very different properties, yet following the general trend that in 



the Z state , QPOs showed lo wer Q and amplitude than in the atoll state (IHoman et al. 



2007. 


2010; 


Sanna et al. 


2010) 



between the Z and atoll state spans a factor of 15-20. The dramatic changes in the QPO 
parameters between the two source states could not be due to changes in the neutron star 
mass, its magnetic field, its spin, or even the inclination of the accretion disk. It is thus more 
plausible that they w e re ins tead caused by a change in the properties of the accretion flow. 



This led ISanna et al.l (120101 ) to conclude first that effects other than the geometry of space 
time around the neutron star have a strong influence on the coherence and amplitude of the 
kHz QPOs, and second that the drop of the coherence and RMS amplitude of the lower kHz 
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QPOs, as we have observed it, could not be used to infer the existence of the innermost 
stable circular orbit around a neutron star. 



In this paper, we revisit the RXTE observations reported by ISanna et al.l ( 120101 ). fo- 
cussing on the lower kHz QPO, and apply th e same analysis procedures (frequency drift 
correction) as described in IB arret et al.l ( 20061 ). The main reason is that, despite very few 
detections overall (12 segments of observations in total), the frequency span of the lower 
kHz QPO ranges from about 640 Hz to 850 Hz, and by comparison with other sources, this 
should be s ufficient to investig ate the drop of its quality factor (n o such drop i s obvio us in 



Figure 3 of ISanna et al.l ( 120101 ) who used a different analysis than iBarret et al. 



(2006)). In 



the ne xt section, we describe the re-analysis of the lower kHz QPOs reported by ISanna et al. 
(120 lOl ) to show that the drop of its quality factor and RMS amplitude is indeed observed. 
We then discuss the changes in the properties of the QPOs from XTE J1701-462 as simply 
related to a likely decrease of the disk thickness between t he Z and atoll phases , as previously 
discussed in the framework of the toy model presented in lBarret et al.l ( 120071 ). 



Observations 



We follow the definition of the Z and atoll states as in lSanna et al.l (120101 ): the separation 
is estimated around the end of April 2007. We have retrieved from the HEASARC archive 
science event mode and single bit data recorded by the RXTE Proportional Counter Array 
(PCA). We consider data as segments of continuous observation (an ObsID may contain 
more than 1 segment). For each segment, we have computed an average Power Density 
Spectrum (PDS) with a 1 Hz resolution and an integration time of 16 s econds, using event s 
recorded between 2 and 40 keV. The PDS are normalized according to iLeahy et al.l ( 119831 ). 
so that the Poisson noise level is expected to be a constant close to 2. The PDS is then 
blindly searche d for excess power b etween 500 Hz and 1400 Hz using a scanning technique, 
as presented in iBoirin et al.l ( 120001 ) . We have also verified that no significant excesses were 
detected between 1400 and 2048 Hz. This justifies the use of the 1400 - 2048 Hz range to 
estimate accurately the Poisson noise level in each observation, which is indeed close to 2 in 
all segments of data. 

The excess power is then fitted with a Lorentzian with three free parameters; frequency, 
full width at half maximum, and amplitude (equal to the integrated power of the Lorentzian). 
The Poisson noise level is fitted separately above 1400 Hz. 68% confidence errors on each 
param eter are computed in a standard way, using a Ax 2 of 1. Following iBoutelier et al. 
(120091 ). our threshold for QPOs is related to the ratio (hereafter R) of the Lorentzian am- 
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plitude to its la erroi0 (R was often quoted and used as a significance). A conservative 
threshold R = 3 means that we consider only QPOs for which we can measure the power 
of the Lorentzian with an accuracy of 3a or more. Such a threshold corresponds to a ~ 6a 
excess power in the PDS for a si ngle trial, equivalen t to ~ 4a if we account for the number of 
trials of the scanning procedure (Ivan der Kliall989h . The integrated power of the Lorentzian 
is then converted into a root mean square (RMS), expressed as a fraction of the total source 
count rate. 



We detect a s i gnific ant QPO in twelve segments of observations, within the same ObsIDs 



as m 



Sanna et al.l ( 120101 ). The reduced number of active PCUs on the PCA, and the relative 
faintness of the source in its atoll phase (count rate between ~ 40 and ~ 120 counts/s/PCU), 
means that special care must be taken to c orrect for the f requen cy drift of the lower kHz QPO. 
We apply the very same technique as in iBarret et al.l ( 120061 ). It is an iterative procedure, 
which enables us to bound the QPO frequencies on shorter and shorter integration time, with 
narrower and narrower frequency interval s. As the integratio n time decreases, the significance 
threshold set to the scanning technique (IBoirin et al.ll2000l ) is adjusted (e.g the threshold is 
decreased from 4<r for an interval of 50 Hz over a few thousand seconds, to 3a for an interval 
of 10 Hz over a few hundreds of seconds). The QPO path is then recovered through a linear 
interpolation between the most significant detections with the shortest integration time. An 
example of the application of the analysis to two segments of real observations is shown in 
Figure CD 

Such a technique, which we have used extensively on real and simulated data, is able 
to correct precisely for the frequency drift of QPOs of the strength reported here. In those 
simulations, the QPO frequency evol ution is modeled by a ran dom walk of a given step 
(e.g. 0.2-0.5 Hz/second). We follow iTimmer &: Koenigl (119951 ) and generate synthetic 1 
second PDS, for which the underlying model consists of a constant (2) to account for the 
Poisson noise plus a lorentzian to model the QPO profile, with parameters (Q and RMS) 
appropriate for each frequency (estimated from the random walk). The simulated PDS 
are then combined on (16 second timescales) and scanned as the real data to recover the 
time evolution of the QPO frequency, and to determine Q and RMS, for comparison with 
the parameters injected in the simulations. The difference between the reconstructed and 



1 The Lorentzian function used in the fit is Lor(z/) = A x W/(2n)/[(i> — vq) 2 + (W/2) 2 ], where A is 
the integrated power of the Lorentzian from to oo, W its width and vq its centroid frequency . The fitted 
function is linear in A, and therefore its error can be computed using A\ 2 (e.g. iPress et al.l (119921)). The RMS 
ampl itude is a derived quantity, computed as RMS=\/A/S, where S is the source count rate ( van der Klisl 
1989 ). In this paper, we have defined R = A/SA, from which the error on the RMS is estimated as 
<5RMS = 1/2 x RMS x R" 1 after neglecting the term SS/S in the derivative of the RMS equation. 
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original Q and RMS parameters can then be evaluated. Those simulations have shown that 
for QPOs of comparable strength to the one of XTE J1701-462 (mean RMS around 9%), 
Q can be recovered with an accuracy better than 10% (the accuracy on the RMS is better 
than a few %). 

Having reconstructed the QPO frequency evolution, one can compute the mean QPO 
parameters over the observations, by aligning all the 16 second PDS for which a QPO 
frequency was estimated, directly from a significant detection or by interpolation between 
two significant detections, separated by less than 256 seconds, (see Figured]). Increasing the 
latter value to say, 512 seconds, does not lead to any changes in the fitted QPO parameters. 
The best fit results are listed in Table [1] and presented in Figure [2j All the QPOs reported 
have a R > 5, and are therefore highly significant. As can be seen from Figure [2] despite 
some scatter, there is already evidence that above 800 Hz the quality factor and the RMS 
amplitude of the QPO drops. 

In addition to those highly significant QPOs, we also note that there are hints for two 
additional single QPOs in the segments 93703-01-03-04 and 93703-01-05-06. No frequency 
drift corrections are possible for those two QPOs. In the first segment, v = 905.0 ± 4.4 Hz, 
Q = 40.6 ± 19.3 for an R ratio of 2.9, while in the second one u = 548.3 ± 33.2, Q = 2.7 ± 1.6 
for an R factor of 2.3. Being single and not very significant, it is difficult to draw any firm 
conclusions. However, it cannot be excluded that they are lower kHz QPOs, extending on 
both sides the frequency span of figure [21 in which case one would expect them to have low 
Q factors, as measured (with indeed the caveat that no drift correction could be applied, 
implying that the measured values are only lower limits on Q). 

The histogram of interpolated frequencies (measured over 16 seconds) is shown in Figure 
El As can be seen, the frequency span of the lower QPO, although comparable in breadth 
with other systems, has not been sampled equally (there is a lack of observations around 
750 Hz). Grouping the data over constant frequency intervals is therefore not the optimum 
way to proceed. Instead we have considered adjacent frequency intervals of varying widths, 
each including 200 QPO frequencies, allowing a better sampling of the peak of the quality 
factor versus frequency curve. The result is shown in Figure SJ As can be seen, the drops 
around 800 Hz of the QPO RMS amplitude and quality factor are now much clearer and have 
significantly less scatter. Thanks to our QPO tra cking proced u re, w e can recover the QPO 



frequency on a timescale of 16 seconds, whereas ISanna et al.l (120101 ) averaged as many 16 



second PDS as required to enable a s i gnific ant detection. This explains why on average we 



get larger Q factors than ISanna et al.l ( 120101 ) . and hence a better description of its frequency 
dependence. 



Following ISanna et al.l (120101 ). we have shifted all the 16 second PDS to search for the 
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upper kHz QPO. ISanna et al.l (120101 ) reported a 3.1<r significance detection (single trial) with 
a frequency separation of 258 ± 13 Hz. Using our procedure, the significance of the upper 
kHz QPO detected rises to 4.8 a (single tria l significance ), for a frequency separation of 
266.5 ± 13.1 Hz, consistent with the value of ISanna et al.l ( 120101 ). At the same time, the 
significance of the main peak increased from 30<r to more than 50er in our analysis (see 
Figure [5]). In the two ObsIDs of the Z phase, in which the simultaneous twin kHz QPOs 
are the most significant (R > 3) (92405-01-40-04 and 92405-01-40-05), we have measured a 
frequency difference of 280.8 ± 17.6 and 276.1 ± 18.8 Hz respectively, i.e. consistent within 
errors with the frequency difference we have measured in the Atoll phase. 



3. Discussion 

Based on the observations reported here, our main findings can be summarized as fol- 



lows: 



• XTE J1701-462 behaves similarly to the other sources we have studied. In the atoll 
state, its lower kHz QPOs are narrow. Correcting for the frequency drift, we have been 
able to show that its quality factor reaches a maximum of 200 around 800 Hz before 
dropping off sharply. This effect is even visible within a limited sample of observations, 
when the drift correction is applied. We have also found a weaker upper kHz QPO 
with a significance greater than previously reported, using the shift-and-add technique. 



Following iBarret et ajj ( 120061 ). we estimate the frequency at the ISCO by adding ~ 270 
Hz 

^lower.c^cb which we infer to be around 900 Hz (albeit with large uncertainty due 
to the limited sample of the Q — v curve in the available data set). This gives i^isco ~ 
1170 Hz. No upper QPO above that f requency should be detected in XTE J1701-462. 
This is consistent with the findings of ISanna et al.l (120101 ). who reported upper QPOs 
with frequencies less than ~ 930 Hz. From z/t sco, we can estimate the gravitational 
masses of the neutron stars (IMiller et al.lll998l ): 



M « 2.2 M o (1000 Hz/z/is C o)(l + 0.75;') 



(1) 



where j = cJ/(GM 2 ) ~ 0.1 — 0.2 is the dimensionless angular momentum of the star. 
The inferred mass of the neutron star in XTE J1701-462 would therefore be ~ 2.0 M . 



• As shown by lSanna et al.l (120101 ). when the source is in the Z phase, the properties of its 
QPOs are consistent with those observed in other Z sources: they are weak and broad 
(the maximum Q values are about 10 and the RMS less than ~ 4%). For the purpose of 
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this paper, we have also re-analyzed the segments of data in which ISanna et all (120 101 ) 
reported a QPO, and our results are globally consistent with theirs. Note however 
that no frequency drift corrections can be applied when recovering the quality factor, 
hence the measured values should be considered as lower limits. In any case, our 
independent analysis confirms that no highly coherent QPOs have been detected from 
XTE J1701-462 in the Z phase. This indeed suggests that there is a para meter which 



alters the Q value, in addition to the radius at which the QPOs originate. ISanna et al. 



(120101 ) conclude that these results weaken the hypothesis relating the drop in coherence 
observed in many sources, and in the Atoll phase of the one under investigation, to the 
presence of the ISCO, because there are evidently other properties of the accretion flow 
that can change the quality factor of QPOs. We point out, though, that the drop in 
coherence that we relate to the ISCO is a very characteristic phenomenon. In a given 
source, it happens always at the same frequency, at the higher end of the frequency 
range of kHz QPOs. There is therefore a significant difference between the drop in 
coherence we relate to the ISCO and the change in coherence in the different states of 
this source. In our opinion, this implies that different phenomena are producing the 
two effects. 



Let us now discuss on some possible interpretations. The change in coherence of the 
kHz QPOs in XTE J1701-462 between the Atoll and Z phase, has to be related to a global 
change of the accretion configuration. 



In the framework of the toy model proposed by lBarret et al.l (120071 ). we argue that the 
difference is related to the scale height of the disk, exactly the same process that we propose 
as responsible for the general difference in coherence between QPOs in Atoll and Z sources. 
Consistent with this idea, standard disk accretion theory suggests that at rates approaching 
Eddington (which is expected in the Z state) the disk half-thickn ess at its inner edge will be 
comparable to the orbital radius there ( iShakura &: Sunyaevlll973l ). and that as a consequence 



the inward radial drift speed (whi ch scales as 
be large as well. As discussed in iBarret et al. 



h/r) 2 , where h is the disk half-thickness) will 
(120061 ) . a large inward speed will necessarily 



decrease Q regardless of other fac tors. It is the refore not surprising that high luminosity 
sources have lower Q, as observed ( IMendezI 120061 ). 



We note that alternative explanations are under study through 3D MHD simulations 
that show oscil l ation s promisingly similar to the kilohertz QPOs (for a recent discussion see 
Bachetti et al.l ( 120101 )). In these simulations, the lower and upper kHz QPOs are emitted from 
hotspots on the surface of the star, produced by streams of matter originating respectively 
from the funnel flow around the magnetic pole and from instabilities at the inner edge of 
the disk. The properties of these QPOs in simulated systems with an enhanced variety of 
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physical inputs (including different configurations of the disk) are currently under study. 



Conclusions 



The source XTE J1701-462, which has been observed in both the atoll and Z phases, 
gives us a unique opportunity to explore the relation between frequency and quality factor. 
This relation has been claimed for previous atoll sources to provide evidence of the innermost 
stable circular orbit. It was expected, however, that for higher-luminosity sources, in which 
the disk becomes geometrically thick and the inward radial speed is thus large, the quality 
factors woul d be systemat ically lower than they are in the lower-luminosity atoll sources, 
as observed (jMendezll2006l ). XTE J1701-462 does show exactly this behavior. In the atoll 
state its Q vs. v phenomenology is consistent with that of other sources, but in the Z state 
the quality factors are systematically smaller than in the atoll state. The change is over 
the whole range of the QPO frequencies and is unlike the abrupt drop observed in the atoll 
phase, which only happens at t he hig her end of the frequency range. As a consequence, we 
do not agree with I S anna et al.l (120 lOh that the different coherence in the atoll and Z phases 
of XTE J 170 1-462 weakens the ISCO hypothesis. 

We can ask: what behavi or would contradict the ISCO picture? As discussed in other 
papers (e.g. iBarret et al.ll2007l ). there are several possibilities including (1) a lower or upper 
kHz QPO frequency in another state that is significantly larger than the projected maximum 
due to the ISCO, (2) a peak in Q followed by a sharp drop at a frequency significantly less 
than seen in the atoll state, or (3) a peak in Q followed by a sharp drop in any source at a 
frequency that would imply an unreasonable ISCO frequency (e.g., below 800 Hz). None of 
these have yet been seen. Therefore, although the importance of the implications demand 
further critical analysis, the interpretation is still consistent with all existing data to date. 

MCM acknowledges NSF grant AST0708424. The authors wish to thank an anonymous 
referee for helpful comments. 
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Fig. 1. — The QPO frequency history recovered by our analysis for two segments of ob- 
servation of the ObsID 93703-01-02-05. The path followed by the QPO is indicated with 
the yellow filled circles. As can be seen on the left, there are gaps in the QPO detection. 
When the gap duration is less than 256 seconds, frequencies within the gap are estimated 
by interpolating between the two positive detections surrounding the gap. 

Table 1. Lower kHz QPOs detected in the atoll phase of XTE J1701-462. Parameters are 

averaged over the segments of observation, after aligning all the PDS to a reference 
frequency. The table lists the segment number, the ObsID name, the date, the duration of 
the segment, the mean count rate normalized by the number of active PCUs, the reference 
frequency of the QPO, its quality factor, its RMS amplitude, and the accuracy R ratio, as 

defined in the text. 



Num 


ObsID 


Date 


Time 


Duration 


C 


V 


Q 


RMS 


R 


01 


93703-01-02-04 


2007/07/24 


08 


11 


11 


1968 


116.9 


849.9 ±0.5 


126.2 ±26.0 


7.8 ±0.6 


7.0 


02 


93703-01-02-04 


2007/07/24 


10 


15 


27 


1280 


117.3 


854.2 ± 0.7 


105.7 ±24.6 


8.7 ±0.7 


6.3 


03 


93703-01-02-05 


2007/07/25 


10 


02 


19 


1920 


105.8 


837.5 ± 0.3 


138.3 ± 24.3 


9.5 ±0.5 


8.8 


01 


93703-01-02-05 


2007/07/25 


11 


23 


27 


3376 


66.7 


832.9 ±0.2 


152.2 ± 13.2 


8.8 ±0.3 


16.0 


05 


93703-01-02-08 


2007/07/25 


13- 


11 


10 


2304 


104.9 


835.4 ±0.2 


176.5 ± 21.8 


9.1 ±0.4 


11.0 


06 


93703-01-02-11 


2007/07/25 


05 


13 


35 


2528 


108.2 


834.9 ± 0.3 


175.1 ± 31.5 


7.7 ±0.5 


8.1 


07 


93703-01-02-11 


2007/07/25 


06 


40 


31 


3392 


68.1 


802.8 ±0.1 


189.9 ± 13.1 


9.4 ±0.2 


20.4 


08 


93703-01-02-11 


2007/07/25 


08 


11 


23 


208 


71.0 


786.7 ± 0.4 


168.7 ±37.2 


10.3 ±0.8 


6.3 


09 


93703-01-03-00 


2007/07/29 


03 


19 


30 


3344 


56.4 


772.7 ±0.8 


94.3 ± 24.8 


10.5 ±0.9 


5.6 


10 


93703-01-03-00 


2007/07/29 


04 


53 


19 


3408 


39.2 


712.1 ± 0.6 


71.3 ±9.4 


11.4 ±0.6 


9.7 


11 


93703-01-03-00 


2007/07/29 


06 


27 


27 


2976 


42.1 


707.7 ±0.4 


103.4 ± 17.1 


10.3 ±0.6 


9.1 


12 


93703-01-03-02 


2007/07/29 


11 


23 


57 


1360 


39.7 


655.1 ± 1.0 


73.6 ± 19.3 


9.9 ±1.0 


5.2 
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Fig. 2. — RMS amplitude (top panel) and quality factor (bottom panel) of the lower kHz 
QPO from XTE J1701-462 as averaged over the twelve segments of observations. There is a 
trend for the quality factor and the RMS amplitude to decrease after reaching a maximum 
around 800 Hz, in a very similar way as observed in other sources. 
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Fig. 3. — Histogram of the frequencies (one frequency per 16 seconds) detected in all 12 
observations, showing that not all frequencies have been sampled equally. The frequencies 
used for building the histogram are available in an online table. 
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Fig. 4. — Same as Figure 2, but the QPO parameters are averaged over frequency intervals 
(200 consecutive frequencies were combined). The drop in Q and RMS amplitudes are now 
clearer. 
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Fig. 5. — The two kHz QPOs detected by shifting and adding all the 16 second PDS in which 
a QPO frequency could be estimated. The main peak has a significance above 50<r, whereas 
the upper kHz QPO has a single trial significance of 4.8 a. Its quality factor is 17.1 ± 9.1 
and its RMS amplitude 4.4 ± 0.8, and its R ratio is 2.6 



